Projeto

Leitura dos dados

dados <- haven::read_sav(here("data", "STU_QQQ_5.sav"))
dados

Omissos de cada coluna

colSums(is.na(dados))

     PV1SCIE      PV1READ      PV1MATH      SOCONPA      BODYIMA     GCAWAREP     INTCULTP 
           0            0            0       179328       180414       199429       203838 
     ATTIMMP     JOYREADP      PRESUPP     PASCHPOL     PQSCHOOL      EMOSUPP      CURSUPP 
      203832       199438       199113       199379       199136       199549       199040 
    FLFAMILY     FLSCHOOL     FLCONICT     FLCONFIN     INFOJOB2     INFOJOB1      INFOCAR 
      141018       139894       149718       148000       136593       136593       130465 
  ICTOUTSIDE     ICTCLASS      SOIAICT       AUTICT      COMPICT       INTICT       USESCH 
      106198       102646       118574       116498       110049       107134       115491 
     HOMESCH       ENTUSE BEINGBULLIED       BELONG      DISCRIM     GLOBMIND       AWACOM 
      107699        97598        55813        27833       116537        84839        62117 
     RESPECT      COGFLEX     PERSPECT      INTCULT       ATTIMM      GCAWARE    GCSELFEFF 
       65712        61036        60271        81288        87447        53427        57851 
    MASTGOAL   RESILIENCE         SWBP        EUDMO      GFOFAIL     WORKMAST      COMPETE 
       27371        30429        46607        49178        27875        33339        22776 
    ATTLNACT      PERCOOP      PERCOMP     PISADIFF   SCREADDIFF   SCREADCOMP      JOYREAD 
       21095        76270        70099        19611        20720        26831        12444 
    TEACHINT   ADAPTIVITY     STIMREAD      EMOSUPS      PERFEED       DIRINS     TEACHSUP 
       34011        15204        14616        66628        35518        32959        31187 
    DISCLIMA       ICTRES       WEALTH       HEDRES     CULTPOSS      HOMEPOS       ICTSCH 
        9475         6996         6319         7687         9285         5887        94372 
     ICTHOME         ESCS       STUBMI       CHANGE     SCCHANGE     FCFMLRTY        TMINS 
       91478         7086       180775       133834       132921       136725        89489 
       SMINS        LMINS        MMINS          AGE        GRADE   ST061Q01NA   ST060Q01NA 
       62120        60320        60475            0          193        50632        61877 
  ST059Q03TA   ST059Q02TA   ST059Q01TA   ST016Q01NA        IMMIG        HISEI        BFMJ2 
       56844        54819        54531        67890        11029        19126        43704 
       BMMJ1       HISCED       FISCED       MISCED       ISCEDL    ST001D01T       REPEAT 
       33669         8637        15197        10183         6869            0         9367 
       PROGN        OCOD3        OCOD2        OCOD1    ST004D01T          CNT 
           0            0            0            0            2            0 

# A coluna dos Paises não tem nulos
sum(is.na(dados$CNT))

[1] 0

dados %>%
  mutate(NumNulos = rowSums(is.na(.))) %>%
  group_by(CNT) %>%
  summarize(TotalNulos = sum(NumNulos),
            NumLinhas = n(),
            PercentNulos = (TotalNulos / (NumLinhas * ncol(dados))) * 100,
            MeanScience = mean(PV1SCIE, na.rm = TRUE),
            MeanRead = mean(PV1READ, na.rm = TRUE),
            MeanMath = mean(PV1MATH, na.rm = TRUE),
            .groups = "drop") %>% 
  arrange(desc(NumLinhas))

Espanha <- dados %>% filter(CNT == "ESP")
Espanha

dim(Espanha)

[1] 35943   104

EspanhaSemiLimpa <- Espanha %>% select(where(~ !all(is.na(.))))
dim(EspanhaSemiLimpa)

[1] 35943    95

EspanhaLimpa <- na.omit(EspanhaSemiLimpa)
dim(EspanhaLimpa)

[1] 8458   95

Espanha <- EspanhaLimpa %>% select_if(is.numeric)
rm(EspanhaLimpa)
rm(EspanhaSemiLimpa)

#Correlation plot with colors (too many attributes)
correlation <- cor(Espanha)
ggcorrplot::ggcorrplot(correlation)

#Correlation matrix
round(correlation, 3)

             PV1SCIE PV1READ PV1MATH SOCONPA BODYIMA FLFAMILY FLSCHOOL FLCONICT FLCONFIN INFOJOB2
PV1SCIE        1.000   0.756   0.740   0.068  -0.048   -0.067   -0.134    0.023    0.022   -0.043
PV1READ        0.756   1.000   0.704   0.089  -0.113   -0.032   -0.056   -0.004   -0.025   -0.053
PV1MATH        0.740   0.704   1.000   0.059  -0.017   -0.077   -0.098    0.017    0.035   -0.040
SOCONPA        0.068   0.089   0.059   1.000   0.182    0.145    0.014    0.010   -0.004   -0.036
BODYIMA       -0.048  -0.113  -0.017   0.182   1.000    0.065    0.038    0.104    0.127    0.035
FLFAMILY      -0.067  -0.032  -0.077   0.145   0.065    1.000    0.188    0.158    0.163    0.137
FLSCHOOL      -0.134  -0.056  -0.098   0.014   0.038    0.188    1.000    0.165    0.205    0.084
FLCONICT       0.023  -0.004   0.017   0.010   0.104    0.158    0.165    1.000    0.744    0.182
FLCONFIN       0.022  -0.025   0.035  -0.004   0.127    0.163    0.205    0.744    1.000    0.198
INFOJOB2      -0.043  -0.053  -0.040  -0.036   0.035    0.137    0.084    0.182    0.198    1.000
INFOJOB1      -0.075  -0.035  -0.051   0.045   0.077    0.093    0.208    0.094    0.105   -0.142
             INFOJOB1 INFOCAR ICTOUTSIDE ICTCLASS SOIAICT AUTICT COMPICT INTICT USESCH HOMESCH
PV1SCIE        -0.075   0.024      0.031    0.106  -0.059  0.100   0.061  0.019 -0.064  -0.111
PV1READ        -0.035   0.068      0.045    0.071  -0.135  0.075   0.062  0.057 -0.086  -0.137
PV1MATH        -0.051   0.037      0.013    0.059  -0.027  0.107   0.064  0.033 -0.056  -0.102
SOCONPA         0.045   0.024      0.046    0.054   0.050  0.017   0.063  0.048  0.004   0.045
BODYIMA         0.077  -0.007      0.002    0.014   0.110  0.072   0.094  0.022  0.056   0.075
FLFAMILY        0.093   0.136      0.134    0.077   0.142  0.109   0.113  0.122  0.127   0.149
FLSCHOOL        0.208   0.150      0.084    0.076   0.090  0.057   0.037  0.012  0.154   0.132
FLCONICT        0.094   0.125      0.067    0.068   0.168  0.212   0.201  0.071  0.117   0.104
FLCONFIN        0.105   0.119      0.064    0.050   0.166  0.183   0.157  0.010  0.101   0.099
INFOJOB2       -0.142   0.180      0.082    0.020   0.137  0.156   0.144  0.038  0.065   0.121
INFOJOB1        1.000   0.201      0.080    0.076   0.104  0.023   0.032  0.018  0.165   0.151
             ENTUSE BEINGBULLIED BELONG DISCRIM GLOBMIND AWACOM RESPECT COGFLEX PERSPECT INTCULT
PV1SCIE      -0.023       -0.059 -0.021  -0.211    0.096  0.208   0.123   0.101    0.063   0.084
PV1READ      -0.054       -0.059 -0.015  -0.262    0.161  0.284   0.207   0.082    0.113   0.182
PV1MATH      -0.033       -0.062  0.012  -0.155    0.093  0.182   0.082   0.071    0.028   0.058
SOCONPA       0.008       -0.191  0.198  -0.169    0.104  0.129   0.125   0.113    0.090   0.063
BODYIMA       0.066       -0.166  0.258  -0.004    0.018  0.007  -0.014   0.187    0.012  -0.052
FLFAMILY      0.142        0.038  0.077   0.053    0.096  0.109   0.037   0.103    0.063   0.121
FLSCHOOL      0.084        0.051  0.048   0.043    0.059  0.047  -0.020   0.067    0.047   0.088
FLCONICT      0.152        0.008  0.053   0.040   -0.007  0.003  -0.040   0.121    0.004   0.019
FLCONFIN      0.126        0.010  0.043   0.051   -0.007 -0.007  -0.058   0.138    0.006   0.002
INFOJOB2      0.141        0.060 -0.004   0.110    0.036  0.007   0.008   0.138    0.078   0.073
INFOJOB1      0.076       -0.012  0.087  -0.010    0.043  0.011   0.008   0.057    0.024   0.050
             ATTIMM GCAWARE GCSELFEFF MASTGOAL RESILIENCE   SWBP  EUDMO GFOFAIL WORKMAST COMPETE
PV1SCIE       0.104   0.061     0.252    0.145      0.093 -0.044 -0.107  -0.003    0.036   0.068
PV1READ       0.173   0.091     0.245    0.164      0.065 -0.067 -0.124   0.060    0.094   0.015
PV1MATH       0.076   0.027     0.205    0.154      0.101 -0.005 -0.087  -0.016    0.055   0.099
SOCONPA       0.069   0.091     0.063    0.220      0.193  0.284  0.204  -0.093    0.217   0.036
BODYIMA      -0.029   0.075     0.094    0.103      0.358  0.326  0.296  -0.313    0.141   0.142
FLFAMILY      0.038   0.167     0.109    0.129      0.124  0.066  0.150   0.041    0.142   0.066
FLSCHOOL      0.009   0.143     0.150    0.070      0.068  0.041  0.116  -0.006    0.065   0.042
FLCONICT     -0.031   0.128     0.167    0.036      0.155  0.055  0.108  -0.113    0.030   0.095
FLCONFIN     -0.070   0.130     0.203    0.040      0.153  0.062  0.124  -0.132    0.035   0.109
INFOJOB2     -0.005   0.148     0.148    0.015      0.099 -0.015  0.075  -0.019    0.044   0.075
INFOJOB1      0.000   0.080     0.059    0.095      0.079  0.089  0.109  -0.032    0.071   0.043
             ATTLNACT PERCOOP PERCOMP PISADIFF SCREADDIFF SCREADCOMP JOYREAD TEACHINT ADAPTIVITY
PV1SCIE         0.055   0.060  -0.017   -0.331     -0.244      0.277   0.240    0.089      0.031
PV1READ         0.086   0.054  -0.033   -0.385     -0.305      0.373   0.372    0.101      0.050
PV1MATH         0.076   0.087  -0.014   -0.288     -0.197      0.206   0.170    0.074      0.034
SOCONPA         0.161   0.186   0.021   -0.063     -0.058      0.054   0.020    0.137      0.132
BODYIMA         0.039   0.157   0.042   -0.039     -0.080      0.031  -0.140    0.058      0.052
FLFAMILY        0.058   0.081   0.115   -0.005     -0.038      0.045   0.035    0.057      0.047
FLSCHOOL        0.010   0.091   0.076    0.011     -0.011      0.066   0.035    0.070      0.063
FLCONICT       -0.033   0.060   0.057   -0.075     -0.072      0.061  -0.053    0.004      0.026
FLCONFIN       -0.029   0.058   0.066   -0.077     -0.082      0.063  -0.059    0.013      0.023
INFOJOB2       -0.022  -0.012   0.081   -0.039     -0.065      0.073   0.021   -0.036     -0.032
INFOJOB1        0.064   0.119   0.046    0.009     -0.014      0.013  -0.020    0.085      0.108
             STIMREAD EMOSUPS PERFEED DIRINS TEACHSUP DISCLIMA ICTRES WEALTH HEDRES CULTPOSS
PV1SCIE         0.012   0.037  -0.050 -0.046    0.015    0.085  0.131  0.072  0.085    0.216
PV1READ         0.010   0.077  -0.068 -0.038    0.035    0.090  0.110  0.033  0.080    0.244
PV1MATH         0.018   0.043  -0.051 -0.048    0.000    0.085  0.167  0.112  0.101    0.210
SOCONPA         0.123   0.531   0.092  0.107    0.149    0.100  0.109  0.115  0.191    0.103
BODYIMA         0.060   0.183   0.104  0.086    0.061    0.060  0.053  0.085  0.098    0.010
FLFAMILY        0.114   0.166   0.119  0.079    0.046   -0.003  0.024  0.032  0.110    0.067
FLSCHOOL        0.128   0.041   0.121  0.072    0.058    0.047  0.000  0.010  0.068    0.036
FLCONICT        0.031   0.024   0.077  0.022   -0.004    0.029  0.040  0.052  0.056    0.033
FLCONFIN        0.052   0.000   0.096  0.030   -0.001    0.029  0.010  0.030  0.051    0.033
INFOJOB2        0.014  -0.017   0.037 -0.016   -0.056   -0.034  0.007  0.013  0.043    0.085
INFOJOB1        0.156   0.062   0.136  0.124    0.106    0.041  0.028  0.026  0.088    0.020
             HOMEPOS ICTSCH ICTHOME   ESCS STUBMI CHANGE SCCHANGE FCFMLRTY  TMINS  SMINS  LMINS
PV1SCIE        0.208 -0.098   0.026  0.265 -0.016 -0.149   -0.106    0.043 -0.003  0.245 -0.012
PV1READ        0.197 -0.111   0.001  0.253 -0.019 -0.142   -0.091    0.115 -0.012  0.165 -0.027
PV1MATH        0.234 -0.083   0.052  0.307 -0.042 -0.171   -0.124    0.065  0.009  0.252 -0.033
SOCONPA        0.165  0.037   0.098  0.145 -0.055 -0.112   -0.101    0.074 -0.012  0.057 -0.010
BODYIMA        0.076  0.083   0.076  0.067 -0.221 -0.033   -0.043    0.058  0.038  0.024  0.029
FLFAMILY       0.067  0.064   0.063  0.021  0.011  0.013    0.000    0.166 -0.007 -0.013  0.001
FLSCHOOL       0.028  0.122   0.026 -0.014  0.010  0.022    0.008    0.470  0.016 -0.265 -0.015
FLCONICT       0.060  0.080   0.068  0.027 -0.001  0.031    0.030    0.220  0.026 -0.025  0.007
FLCONFIN       0.038  0.079   0.045 -0.011  0.026  0.026    0.019    0.252  0.019 -0.033  0.021
INFOJOB2       0.052  0.043   0.047  0.003  0.036  0.083    0.073    0.121 -0.010 -0.011  0.022
INFOJOB1       0.038  0.155   0.064  0.003  0.002 -0.022   -0.037    0.146  0.045 -0.047 -0.020
              MMINS    AGE  GRADE ST061Q01NA ST060Q01NA ST059Q03TA ST059Q02TA ST059Q01TA
PV1SCIE      -0.020  0.058  0.359      0.005     -0.012      0.247     -0.032     -0.018
PV1READ      -0.025  0.060  0.341      0.012     -0.023      0.165     -0.042     -0.037
PV1MATH      -0.033  0.055  0.392      0.008      0.001      0.255     -0.049     -0.042
SOCONPA      -0.019  0.003  0.110     -0.010     -0.007      0.058     -0.012     -0.003
BODYIMA       0.013  0.024 -0.016      0.009      0.030      0.024      0.011      0.028
FLFAMILY     -0.002  0.023  0.004     -0.011     -0.002     -0.015      0.005      0.005
FLSCHOOL     -0.001  0.005  0.032      0.000      0.020     -0.268      0.000     -0.016
FLCONICT      0.029  0.012 -0.015      0.063     -0.003     -0.033     -0.005     -0.019
FLCONFIN      0.020  0.001 -0.032      0.044     -0.002     -0.039     -0.003      0.002
INFOJOB2      0.008 -0.001 -0.068     -0.005     -0.008     -0.011      0.015      0.027
INFOJOB1      0.004 -0.012  0.076      0.061      0.028     -0.053     -0.035     -0.046
             ST016Q01NA  IMMIG  HISEI  BFMJ2  BMMJ1 HISCED FISCED MISCED ISCEDL ST001D01T REPEAT
PV1SCIE           0.033 -0.088  0.238  0.196  0.217  0.173  0.141  0.170  0.032     0.359 -0.381
PV1READ          -0.005 -0.080  0.240  0.205  0.218  0.155  0.131  0.155  0.039     0.341 -0.360
PV1MATH           0.058 -0.109  0.279  0.238  0.258  0.203  0.166  0.206  0.028     0.392 -0.410
SOCONPA           0.325 -0.109  0.090  0.076  0.077  0.084  0.087  0.083  0.008     0.110 -0.125
BODYIMA           0.392 -0.052  0.046  0.039  0.032  0.039  0.047  0.045 -0.018    -0.016  0.015
FLFAMILY          0.073 -0.011 -0.019 -0.020 -0.027  0.007  0.001 -0.001 -0.012     0.004 -0.003
FLSCHOOL          0.042  0.008 -0.043 -0.029 -0.037 -0.012 -0.004 -0.014 -0.013     0.032 -0.036
FLCONICT          0.074  0.053  0.000  0.001 -0.007  0.011  0.013 -0.004 -0.012    -0.015  0.014
FLCONFIN          0.070  0.048 -0.040 -0.029 -0.045 -0.018 -0.020 -0.038 -0.006    -0.032  0.026
INFOJOB2         -0.023  0.067 -0.039 -0.030 -0.035 -0.003  0.006 -0.010 -0.002    -0.068  0.061
INFOJOB1          0.113 -0.015 -0.021 -0.017 -0.009 -0.002  0.001 -0.004  0.010     0.076 -0.078
             ST004D01T
PV1SCIE          0.056
PV1READ         -0.139
PV1MATH          0.078
SOCONPA         -0.074
BODYIMA          0.189
FLFAMILY        -0.053
FLSCHOOL        -0.010
FLCONICT         0.132
FLCONFIN         0.182
INFOJOB2         0.076
INFOJOB1         0.048
 [ reached getOption("max.print") -- omitted 79 rows ]

#KMO measure 
EFAtools::KMO(correlation)

Error in EFAtools::KMO(correlation) : 
  ⓧ Correlation matrix is singular, no further analyses are performed.